Journal of the Physical Society of Japan 



LETTERS 



Dynamical Mean-Field Study of Local Pairing Interaction Mediated by 
Spin and Orbital Fluctuations in Iron Pnictide Superconductors 
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We investigate the two-orbital Hubbard model, which reproduces the electron and hole Fermi surfaces in 
the iron pnictide superconductors, in the presence of the Jahn- Teller electron-phonon coupling by using the 
dynamical mean-field theory. When the intra- and inter-orbital Coulomb interactions, U and U', increase 
with U — U', both the local spin and orbital susceptibilities, Xs an d Xo, increase with Xs = Xo because 
of the spin-orbital symmetry. Due to the Hund's rule coupling J, Xs is enhanced and dominates over Xo 
resulting in the repulsive local pairing interaction VJ oc > 0, while due to the electron-phonon coupling g, 
Xo is enhanced and dominates over Xs resulting in the attractive one V\ oc < which induces the intra- 
orbital s-wave pairing. Remarkably, V\ oc is weakly dependent on doping and can be attractive for heavily 
electron-doped regime where the superconductivity is observed without Fermi surface nesting. 

KEYWORDS: iron pnictide superconductor, multi-orbital system, Jahn-Teller electron-phonon coupling, orbital 
fluctuation, magnetic fluctuation, dynamical mean-field theory 



The iron pnictide superconductors exhibit the com- 
mon feature of phase diagrams, where parent compounds 
show the tetragonal-orthorhombic structural transition 
and the stripe-type antiferromagnetic (AFM) transition 
both of which are suppressed by carrier doping x result- 
ing in the high-T c superconductivity. 1,2 ) When approach- 
ing the AFM transition, the AFM fluctuation observed 
by the NMR experiments 3 ) is found to be enhanced, 
while, when approaching the structural transition, the 
dxz-dyz ferro-orbital (FO) fluctuation (or the O x 2_ y 2 fer- 
roquadrupole fluctuation) 4 ) responsible for the softening 
of the elastic constant Cqq observed by the ultrasonic 
experiments 5, 6 ) is found to be enhanced. Correspond- 
ingly, two distinct s-wave pairings: the s±-wave with sign 
change of the order parameter between the hole and the 
electron Fermi surfaces (FSs) mediated by the AFM fluc- 
tuation 7, 8 ) and the s++-wave without the sign change 
mediated by the FO fluctuation 9, 10 ) and by the antiferro- 
orbital (AFO) fluctuation 11 ) which is also responsible for 
the softening of Cqq through the two-orbiton process, 12 - 1 
were proposed. 

Recent experiments have revealed that the high-T c su- 
perconductivity is realized even in the case with heav- 
ily electron-doped compounds such as RFeAsOi-^Hj; 
(R=Sm, Ce, La) 13 " 15 ) up to x ~ 0.5 and k x Ye 2 ^2 
(A=K, Cs, Rb) 16 ' 17 ) where the large electron FSs are 
observed without the hole FSs. In this case, the mech- 
anisms based on the AFM 7,8 ) and the AFO 11 ) fluctua- 
tions, which are enhanced due to the nesting between the 
electron and hole FSs, seem to be insufficient for explain- 
ing the superconductivity. As for the mechanism based 
on the FO fluctuation 9 ) which is enhanced due to the cou- 
pling between the d xz -d yz orbital fluctuation and the or- 
thorhombic phonon, the superconductivity does not need 
the FS nesting effect but is restricted near the tetragonal- 
orthorhombic structural transition with small x within 



the random phase approximation (RPA) . 

Generally, in the magnetic fluctuation mechanism, the 
the pairing interaction V(q) with wavevector q is repul- 
sive and then the strong q-dependence of V(q) realized 
near the magnetic ordered phase is crucial for the super- 
conductivity. On the other hand, in the orbital fluctua- 
tion mechanism, V{q) is attractive and then the strong 
q-dependence of V(q) realized near the orbital ordered 
phase is not necessary for the superconductivity. When 
the local component of the orbital fluctuation is rela- 
tively larger than that of the magnetic fluctuation, the 
local component of the pairing interaction Vj oc , which is 
nothing but the q-averaged value of V(q), becomes at- 
tractive and induces the s ++ -wave pairing, even far away 
from the ordered phases. 

In this letter, we investigate the local paring interac- 
tion V\ oc mediated by the spin and orbital fluctuations 
by using the dynamical mean- field theory (DMFT) 18 ) 
which becomes exact in infinite dimensions {d — oo) and 
enables us to sufficiently include the local correlation ef- 
fect beyond the RPA. We employ the two-orbital Hub- 
bard model, 19,20 ) which reproduces the electron and hole 
FSs in the iron pnictides, in the presence of the coupling 
between the orbital fluctuation and the Jahn-Teller (JT) 
phonon corresponding to the orthorhombic mode respon- 
sible for the softening of Cqq as discussed in the previous 
work. 9, 10 ) Although the present model is a simplified ver- 
sion of the multi-orbital electron-phonon models for the 
iron pnictides, 9 ~ 12 > 21 ) the essential feature of the local 
magnetic and orbital fluctuations, which is crucial for 
determining V\ oc , is expected to be well described. 

Our model Hamiltonian is given by 



with the kinetic part of the Hamiltonian: 

dkla 



(1) 



Zfc2cr 



H k =( 6kl ekl2 ) (2) 
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where d^ia is the annihilation operator for a Fe 3d elec- 
tron with the wavevector k and the spin a in the orbital 
I = 1, 2(= d xz ,d yz ), and the energies tki and £ki2 are 
determined so as to reproduce the electron and hole FSs 
in the iron pnictides. 19 ^ The Coulomb interaction part 
flint includes the intra- and inter-orbital direct terms U 
and U' , the Hund's rule coupling J and the pair transfer 
J'. For simplicity, we assume the relations U = U' + 2J 
and J = J' which are satisfied in the isolated atom but 
not generally in the crystal. 22,23 - 1 The phonon and the 
electron-phonon interaction parts are given by 



(3) 



where bi is the annihilation operator for a JT phonon at 
site i with the frequency uq, which is coupled to the lon- 
gitudinal orbital fluctuation, r Z i = 'Yl l(7 {nii a — n^o-) with 
nua = dlia-dua, through the electron-phonon coupling g. 

To solve the model eq.(l), we use the DMFT 18 ) in 
which the lattice model is mapped onto an impurity An- 
derson model embedded in an effective medium which is 
determined so as to satisfy the self-consistency condition 



g-\ z )-E(z) 



z — Hk- 



(4) 



where S(z) and 0(z) are the 2x2 matrix representations 
of the impurity (local) self-energy and the bare impurity 
Green's function, respectively, and z is the complex fre- 
quency. 

We write the effective pairing interaction for the spin- 
singlet state mediated by the spin and charge-orbital fluc- 
tuations in the 4x4 matrix representation as 24 ) 



V(q,z) = -T s (z)xs(q,z)T s (z) 



1 



i ( f(°) + rf)), 



-T c (z)xc{q,z)T c (z) 
(5) 



where T s m(z) is the local irreducible spin (charge- 
orbital) vertex in which only the external frequency 
(z) dependence is considered as a simplified approxi- 
mation. 25 ) The corresponding bare vertex is given by: 
F$)]hh = U(U-2g 2 D(z)), [T^ c) } wl V = U'(—U' + 2 J), 

F&HJ'J' = J ( 2U ' - J + 2 9 2 D{z)) and [T^ c) ]unn 



J'(J') with the bare phonon Green's function D(z) — 
2wo/(w 2 — z 2 ), where V ^ I and the other ma- 
trix elements are 0. 22,23) Within the DMFT, the 
spin (charge-orbital) susceptibility can be expressed 
as Xs(c){q,z) = [Xoiq^y 1 - (+)f s(c) (z)]- 1 , where 
Xo(q, iu) = —T J2kt G(k + q,ie + iuj)G(k, ie) with the 
lattice Green's function G(k, z) — [z — — E(z)]^ 1 , and 
the local vertex can be expressed as T s ^(z) = Xo 1 ( z ) ~ 
Xs( c )( z )' where Xo(^) = -TJ2 £ G(ie + iu;)G(ie) with 
the local (impurity) Green's function G(z) = \Q^ 1 (z) — 
E(z)] -1 . When we replace r s i c )(z) with and ne- 
glect t(z), eq. (5) yields the RPA result of V(q, z). 21 ~ 24 ) 
Therefore, eq.(5) is a straightforward extension of the 
RPA result of V(q, z) to include the vertex and the self- 
energy corresctions within the DMFT without double 
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Fig. 1. (Color online) U(= U') dependence of Z, (S 2 ), (8n 2 ), \s 
and \o with J = J' = g = for n = 2 (lines) and n = 2.2 (marks). 



counting. 

The superconductivity can be examined by solving the 
Eliashberg equation with V(q,z) given in eq.(5) where 
the 1st and 2nd terms of r.h.s. yield fc-dependent anoma- 
lous self-energies (or gap functions) which includes 1/d 
corrections for the DMFT. 18 ) In this letter, we focus 
only on the intra-orbital part of the local paring interac- 
tion V loc = [V(z = 0)]mi with V(z) = ^E q V(q,z) = 

§f s (z)Uz)t s (z) - if c (z)Uz)t c (z) + i(f i 0) +f i 0) ),as 
it is the most dominant contribution for the s++-wave 
pairing due to the orbital fluctuation. 9, 10 ) 

In the actual calculations with the DMFT, we solve 
the effective impurity Anderson model coupled to the JT 
phonon by using the exact diagonalization (ED) method 
for a finite-size cluster to obtain the local quantities 
such as £(z) and Xs(c){ z ) at T = 0, where we set the 
site number N s =4 — 6 and the cutoff of the phonon 
number Nf, = 20. 26,27 ' The tight-binding parameters 
are set to be the same in Ref. 19 ) where the total band 
width is W = 12, and we set the phonon frequency 
ojo = 0.01W. Using the ED method, we calculate the 
several physical quatities as follows: the renormaliza- 
tion factor corresponding to the inverse effective mass 
Z = (1 — d |^ =Q )~ 1 = (m* I m) _1 , the local charge, spin 
and orbital fluctuations (Sn 2 ) with Sn — n — (n), (S 2 ) 
with S = | Y,i J2 a /3 Aia^apduij and (r 2 ), and the local 
spin and orbital susceptibilities Xs = ^((S z \S z ))\ z =o and 
Xo = ((t z \t z ))\z=o, together with the intra-orbital local 
paring interaction V\ oc mentioned above. 

Fig. 1 shows the several physical quantities mentioned 
above as functions of U(= U') with J = J' = g = 
at half- filling n — 2 and away from half- filling n = 2.2. 
When the electron correlation increases with U = U', 
Z and (8n 2 ) decrease while (S 2 ) and (t 2 ) increase with 
(S 2 ) = |(r 2 ) as the double- occupancy probabilities take 
the same value: {ni-\ni\) = {ni^nv\) = (riit n ''t) wit h ' 
V because of the spin-orbital symmetry. 27 ' Correspond- 
ingly, Xs and Xo increase with Xs = Xo while the charge 
susceptibility decreases (not shown) with increasing U. 
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Fig. 2. (Color online) 2g 2 /uiq dependence of Z (a), X s u,) 3 ) an< i 
Hoc ( c ) f° r several values of J(= J') with {/' = 4 and n = 2. 



Fig. 3. (Color online) J dependence of Z (a), \ s and \ (b) and 
Hoc ( c ) f° r several values of 2g 2 /ujq with U' = 4 and n = 2. 



For 7i = 2, we also observe the Mott metal- insulator 
transition at a critical interaction U c (= U' c ) ~ 2.5W, 
where Z = for U > U c , while when U — > U c for 
U < U c , Z — > and x s = Xo — > oo (not shown), as pre- 
viouly observed in the multi-orbital Hubbard model. 28 ' 
For n = 2.2, the U dependence of the physical quantities 
is almost similar to that for n — 2 as shown in Fig. 1, ex- 
cept for the Mott transition which is observed exclusively 
for integer fillings. 28 ' We note that, althogh the q depen- 
dence of X s {q, z) largely dependes on doping responsible 
for the Fermi surface nesting, the q-averaged value, i. 
e., the local susceptibility Xs{z) is weakly dependent on 
doping. 

Next, we consider the effects of the Hund's rule cou- 
pling J(= J') and the JT electron-phonon coupling g. In 
Figs. 2 (a)-(c), Z, x7(o) an< ^ ^ oc are Plotted as functions 
of 2g 2 /ujo for several values of J with U' = 4 and n = 2. 
When 2g 2 /loq increases, Z decreases with increasing Xo 
due to the strong orbital-lattice coupling effect, while Xs 
slightly decreases. Correspondingly, V\ oc decreases with 
increasing 2g 2 /loq and finally becomes negative where the 
attractive term due to Xo dominates over the repulsive 
term due to Xs (see eq.(5)). Then, the intra-orbital s- 
wave (s_)-+-wave) pairing is expected to be realized in 
the intermediate coupling regime where V\ oc < and 
Z = 1/2 ~ 1/3. In the strong coupling regime, we also 
observe the bipolaronic transition at a critical coupling 
g c where Z — > with g — > g c together with x Q — > oo 
and Vi oc —¥ oo, although it is difficult to obtain a fully 
converffent solution with Z ^ 0. The effect of J enhances 



Xs while suppresses Xo (see also Fig. 3). 

Figs. 3 (a)-(c) show the J dependence of Z, xjr ) 
and Vioc for several values of 2g 2 /uJo with U' = 4 and 
n = 2. When J increases, Z monotonically decreases for 
2g 2 /u>o = while it shows a maximum at J ~ 2g 2 /gjq for 
2g 2 Iujq 7^ 0. Xs{Xo) increases (decreases) with increas- 
ing J resulting in a crossover between the following two 
regimes: J S 2g 2 /ujq with Xs < Xo and J Z 2g 2 /uiq with 
Xs > Xo- Then, the large effective mass m* jra = Z~ x 3> 
1 is observed in the two distinct regimes with J <C 2g 2 /uj a 
(J 2g 2 /u>o) where Xo(Xs) dominates over Xs(Xo) due 
to the strong coupling (correlation) effect, while the mod- 
erate effective mass m*/m = Z^ 1 = 2 ~ 3 is observed 
in the intermediate regime with J ~ 2g 2 /u>o where Xs 
and Xo, both of which are largely enhanced by U(= U') 
as shown in Fig. 1, compete to each other resulting in a 
maximum of Z as a fully non-perturbative effect. This 
intermediate regime with Xo being a little larger than 
Xs is responsible for the s++-wave pairing with V\ oc < 
and seems to be consistent with the iron-pnictide super- 
conductors where both the spin and orbital fluctuations 
are large while the renormalization of the band width is 
moderate Z = 1/2 - 1/3. 2 ' 

Finally, we study the doping dependence of V\ oc . In 
Fig. 4, we plot T4oc as a function of n for several values 
of 2g 2 /ujQ with U' = 4 and J = 0.4. For small 2g 2 /uj , 
Vioc is repulsive as Xs dominates over Xo- When 2g 2 /u>o 
increases, Xo increases and then dominates over Xs result- 
ing in V\ oc < as shown in Fig. 2 (c) for n = 2. Remark- 
ably, Vjoc is weakly dependent on n and can be attractive 
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Fig. 4. (Color online) n dependence of Vi oc for several values of 
2g 2 /u with V = 4 and J = J' = 0.4. 



even far away from half-filling, where the s++-wave su- 
perconductivity is expected to be realized without the 
nesting between the electron and hole FSs as recently 
observed in the heavily electron-doped compounds. 13 ~ 17 ) 

In summary, we have investigated the local pairing in- 
teraction Vi oc mediated by the spin and orbital fluctu- 
ations in the two-orbital Hubbard model for iron pnic- 
tides by using the DMFT+ED method which enables us 
to sufficiently include the local correlation effects due to 
the Coulomb and JT electron-phonon interactions be- 
yond the RPA. It has been found that V\ oc is weakly 
dependent on doping and becomes attractive in the inter- 
mediate regime with J < 2g 2 /ujq where the local orbital 
susceptibility is a little larger than the local spin suscep- 
tibility and the band width renormalization is moderate 
1/2 ~ 1/3. The result is responsible for the s++-wave 
pairing even far away from the half-filling without Fermi 
surface nesting and seems to be consistent with the high- 
ly superconductivity observed in the heavily electron- 
doped compounds. 13 ~ 17 ) 

Based on the present study, we can calculate the q- 
dependent spin (charge-orbital) susceptibility X s ( c )(q,z) 
and determine the instability towards the ordered phase 
when the corresponding susceptibility diverges. Prelim- 
inary results shows that the stripe-type AFM observed 
for J > 2g 2 /ujq is largely suppressed as compared to the 
RPA result due to the strong correlation effect within 
the DMFT while the suppression of the FO observed for 
J $ 2g 2 /uj is not so large. We also obtain some re- 
sults for the superconductivity by solving the Eliashberg 
equation and find that the s±-wave occurs near the AFM 
while the s++-waves occur near the FO, where the s++- 
region is qualitatively consistent with the Vi oc < region 
shown in this letter but expands to a rather smaller value 
of 2g 2 /wo. Detailed calculations for the present model, 
together with those for a more realistic 5-orbital model, 
are now under way. 
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